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Abstract. Wc provide a refined approach to the classical Magnus |Mag54| and For expansion IF58I . un- 
veiling a new structure by using the language of dendriform and pre-Lie algebras. The recursive formula for 
. the logarithm of the solutions of the equations X = 1 + Xa ^ X and Y = 1 — XY y a in A[[X]] is provided, 

' where (A,^,^) is a dendriform algebra. Then, we present the solutions to these equations as an infinite 

1 product expansion of exponentials. Both formulae involve the pre-Lie product naturally associated with the 

fS) ' dendriform structure. Several applications are presented 

a ' 

1—5 ; 

Contents 

1. Introduction 
' 2. Dendriform power sums expansions 

■ 3. The pre-Lie Magnus expansion 

f~| . 4. The pre-Lie For expansion 

, 5. AppUcations 

5.1. Associative algebras 

5.2. Rota-Baxter algebras 
Cf^ ■ 5.3. The weight zero case and the classical Magnus expansion 

. References 

o 
p 

. 1. Introduction 

O , 

I Let us start by emphasizing that the results presented in the sections following this introduction arc an 

^ ' extension of findings obtained by the authors together with Frederic Patras in an earher work [EMP07bj . 

The underlying theme of our paper is to bring together Magnus', Fer's and Baxter's classical work on 
, linear differential equations respectively the corresponding integral equations, using the language of Loday's 
dendriform algebras. Skipping details, in this introduction we try to sketch briefly the general picture behind 
our results. 

Recall that Magnus |Mag54| and Fer [F58j , as well as Baxter [Bax60| , start their papers by recalling the 
classical initial value problem: 

(1) ^t):^^^m = ^it)m, m = i- 

Magnus considers it in a non-commutative context, i.e. A .= and Y := <& are supposed to be linear 
operators depending on a real variable t. Here, 1 denotes the identity operator. For the linear operator Q{t) 
depending on A and with f7(0) = 0, such that: 
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y(t) = exp / n{s)ds 







^ nl 

n>0 
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Magnus obtains a differential equation leading to the recursively defined expansion named after him: 

JQ JQ JCl 

Here, as usual, adf[g] := fg — gf [/, ff]- The coefficients i?„ are the Bernoulli numbers defined via the 
generating series: 

1. = y ^ - + _ J_^4 ^ . . . ^ 

exp z - 1 ^ m! 2 12 720 

' m>0 

Observe that i?2m+3 = 0, m > and recall for later use: 

exp(z) — 1 



/ exp(sz) ds. 
Jo 



For more details see for instance |Mag54[ [GKLLRT951 HMNZOOl IK0891 IMFfOl [QROOl [St87l IWUGT] . Observe 
also that Magnus' expansion ([2]) reduces to J* A{s) ds if all commutators disappear, e.g. in a commutative 
setting, leading to the classical exponential solution: 



(3) 



Y(t) = cxp(^ A{s)ds 



for the initial value problem ([T]). Indeed, the reason for this exponential solution is simply encoded in the 
integration by parts rule for I{A){t) := A{s) ds: 



{l{A){t))"=n\l[Al{A---I{A)---)){t). 



n-times 

Fer's |F58j approach to solve the classical initial value problem ([T|), which was rediscovered by Iserles |Ise84| 
and further explored by Zanna [Z96j . and Munthc-Kaas and Zanna [MZ97| . see |IN99j for more details, is 
somewhat different. His Ansatz is: 

Y = exp^ A{s)ds^V{t), 
which leads to the following differential equation for the -correction- operator V{t): 

k>0 

Hence, an iteration of this method leads to the Fer expansion: 

(4) y(t)=exp(^^ Uo{s) dsj exp(^J^ U[(s) ds^ exp(^j^ [/^(s) ds) • • • exp(^ [/^(s) ds) • ■ • , 

where we use Fer's recursion: 

k>0 

Taking an algebro-combinatorial perspective on these methods, we should underline at this point that 
we completely skip the analytical and numerical aspects of these expansions, which are beyond doubt of 
crucial importance in applications. For this purpose and related aspects we refer the interested reader to the 
aforementioned references, e.g. see [BCOR981 IIMNZO"Ol lOROO] . However, at the end of this work we report 
on an observation which may be of interest in this context and which we plan to further explore in the near 
future. 

Baxter |Bax60j considers ([1]) in a commutative setting, i.e. for continuous scalar functions a := and 
y := ^ depending on t. However, his starting point is the corresponding integral equation: 

(6) y{t) = 1+ / a{s)y{s)ds, 

Jo 



(5) t^™+iW-(E7fcTW"4^^^ U[,{t) := A{t). 
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and its exponential solution. Slightly deviating from Baxter's original approach we generalize ([6]) to a formal 
power series ring, VF[[A]], where is a commutative unital algebra over a field k with a fc-linear map 
R : W W replacing the integral map: 

(7) Y =1 + XR{aY), 

a & W fixed. Here, 1 is the unit in the algebra W. The map R is supposed to satisfy the relation: 

(8) Rix)R{y) = R{R{x)y + xR{y) + Oxy) , 

where the parameter is a fixed scalar in fc, called the weight of R. One may think of ([8]) as a generalized 
integration by parts identity. Indeed, the reader will have no difficulty in checking duality of ([5]) with the 
'skewderivation' rule: 

d{fg)^d{f)g + }d{g) + ed{f)d{g). 

For example, the finite difference operator of step —6, given by df{x) :— 6^^{f{x~9) — f{x)), is a skewderiva- 
tion. On a suitable class of functions, the summation operator: 



(9) S{f){x):^Y.^f{x + en). 

satisfies relation (151). Moreover: 



n>l 



Sd{f) {x) = ^ed{f){x + On) ^2^0 

n>l n>l 

n>l ?i>0 Ti>l 

And as the operator d is linear we find as well dS{f) = f. Obviously, the skewderivation rule reduces to 
Leibniz' rule for 6 = 0, e.g. see |CEFG06] for more details. 

The exponential solution of equation jT]) in M^[[A]] can be seen as a natural generalization of the classical 
exponential solution ^ taking into account the weighted term in identity ([5]): 

(10) r ^ i + ^A" i?(ai?(a...i?(a))...) =cxp(i?(i^iii±^)). 

n>0 V ' 

n~times 

The second equality is generally known as Spitzcr's classical identity, see |At63[ lBax601 |EGM06[ IRot69| . In 
fact, expanding the logarithm and the exponential on the right hand side, it follows by comparing order by 
order in the parameter A the infinite set of identities in W^[[A]]. 

Remember that the Riemann integral map / := satisfies identity ([H]) for the weight 6 = (integration 
by parts). In this particular case, observe that: 

6-^ log(l + 6aX) = - V \ a\r ^ aX. 

n>0 

Hence, the exponential on the righthand side of identity (fTO| reduces to the classical exponential solution ([3]). 

Having Magnus' work in mind it seems natural to ask for a non-commutative version of Spitzer's classical 
identity. In fact, pop is only true when the underlying algebra is commutative, i.e. when ([5|) implies: 

R{af = 2R{aR{a)) + dRia^). 

Its generalization to arbitrary weight 6 Rota-Baxter algebras, i.e. non-commutative algebras with a map R 
satisfying relation ([8]) can be found in |EGM06j . Similarly to Magnus' expansion ([2]), it relies on a particular 
non-linear recursion xe, which we call weight 6 BCH-recursion as it is based on the Baker-Campbell- 
Hausdorff formula: 

(11) xe{a)=a+]-BC}i{6a,Roxe{a)). 
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We then find the non-commutative Spitzer identity: 



(12) 




log(l + OaX) 




) 




n>0 « 

n— times 



We refer the reader to Section [5] for more details. One may also consult [EGM061 ICEFGOB] . 

At first sight, the 6 BCH- and Magnus recursion, (fTT|) respectively ([2]), look different, but it is the goal 
of this work to show how they are related. Indeed, in |EGM06) it was already shown that (fTT|) reduces to 
Magnus' formula ^ in the corresponding limit 9^0. Here, on the contrary, we will show how to get pT|) 
from the limit case = 0. After this has been achieved we provide a refined picture of Fer's expansion ([5]) 
in purely algebraic terms leading to an infinite product expansion of the solution to ([7]). These results are 
achieved in the context of the generalized integration by parts rule ([8|) . We emphasize the use of Loday's 
dendriform algebras |Lod01| which appears to be well-suited. 

Recall that a dendriform algebra is an -associative- algebra with two non-associative operations, written 
-< and >- satisfying three rules. The two products add to form the product of the algebra. At the same time 
they define a left and right pre-Lie product on the same algebra. J.-L. Loday recently introduced this notion 
in connection with dialgebra structures [Lod01|, ILR04] . The commutative version of dendriform algebra is 
called "dual Leibniz algebra" or "Zinbiel algebra" by Loday |Lod95| (see also jSchSSp . Temporarily, this 
algebra respectively its product had been called "chronological algebra" respectively "chronological product" 
in the context of control theory, e.g. see [KawOO] for more details. The key point from our perspective is 
the intimate relation between associative algebras, equipped with a map satisfying relation ([8]) and such 
dendriform algebras. We will see that this connection renders dendriform algebras a suitable setting to 
encode algebraic structures related to the integral equations corresponding to linear differential equations. 
At the end we will indicate that this approach unveils a hitherto hidden structure in Fer's and Magnus' original 
work leading to a reduction in number of commutator terms. The reader may also look into [CEFG06] for 
related issues. 

This work is organized as follows. In Section [2] we recall the notion of dendriform algebra and introduce 
two particular dendriform power sums expansions. Section [3] contains the first main result of this work, that 
is, we present the central object, a new pre-Lie Magnus type recursion, and show that its exponentiation 
solves the aforementioned pair of dendriform power sums expansions. Then, in Section [4l we introduce the 
pre-Lie Fer recursion as the second main result, leading to an infinite product expansion of exponentials for 
these solutions. Finally, in Section [5] we provide some applications. We finish the article with an observation 
indicating hitherto overlooked new structural properties of the classical Magnus and Fer expansions due to 
the extra pre-Lie relation, which leads to a reduction in the number of terms in these expansions. 



Let be a field of characteristic zero. Recall that a dendriform algebra |Lod01| over fc is a fc-vector space 
A endowed with two bilinear operations -< and >- subject to the three axioms below: 



In the commutative case, the left and right operations are further required to identify, so that x )^ y — y < x. 
One can show that these relations yield associativity for the product 



2. Dendriform power sums expansions 



(13) 
(14) 
(15) 



(a ^ 6) -< c 
a ^ (6 >- c) 



a ~< (b* c) 
(b ^ c) 
{a*b) >- c. 



(16) 



a * b := a < b + a b. 
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Example 1. As a guiding example we regard an algebra F of operator- valued functions on the real line, 
closed under integrals , say, smooth nxn matrix-valued functions. Then, Dp — {F, -<, >-) is a dendriform 
algebra for the operations: 



{A < B){x) A{x) ■ J B{y) dy {A ^ B){x) J Aiy) dy • 



with A, B £ F. One verifies the dendriform axioms using the integration by parts rule. For instance, in this 
setting the dendriform relation psp means: 

XX X u u 

[a{x) ■ J B{u) du^ ■ J C{v) dv = A{x) ■ J (^B{u) ■ J C{v) dv + J B{v) dv ■ C(it) j du 



for A, B,C £ F. The associative product in this dendriform algebra then writes: 

X X 

{A * B){x) := A{x) ■ I B{y) dy + f Aiy) dy ■ B{x). 



Let us remark that a commutative algebra {F, J^) naturally provides a commutative dendriform algebra. 

Let us return to the dendriform axioms. One shows that at the same time the dendriform relations imply 
that the bilinear operations [> and < defined by: 

(17) a> b -.^ a >~ b — b -< a, a <\b := a ^ b — b >- a 
are left prc-Lic and right prc-Lic, respectively, which means that we have: 

(18) {a> b) \> c- a[> {b\> c) = (6 [> a) [> c - 6 [> (a [> c), 

(19) {a<\b) <lc- a<\{b<lc) = {a < c) <\ b ~ a < {c <i b). 

In the setting of our guiding Example [T] these pre-Lie products write: 

X X 

{A < B){x) := A{x) ■ J B{y) dy - J B{y) dy ■ A{x), 





{A>B){x) :- j A(y)dyB{x) - B{x) ■ J A{y)dy. 



The associative operation * and the pre-Lie operations O, <l all define the same Lie bracket: 

(20) [a, 6] := a *b — b * a = a[> b — b[> a = a <ib — b <\ a. 

Loday and Ronco introduced in [LR04| the notion of tridcndriform algebra T equipped with three oper- 
ations. <, > and •, satisfying seven dendriform type axioms: 

(21) (x < y) < z — X < {y * z), {x > y) < z = x > {y < z), {x -k y) > z ~ x > (y > z), 
(x > y) • z = X > {y • z), {x < y) • z = x • {y > z), {x • y) < z ^ x • {y < z), [x • y) • z ~ x • {y • z), 

yielding an associative product x^y := x < y + x > y + x»y. First, observe that the category of dendriform 
algebras can be identified with the subcategory of objects in the category of tridcndriform algebras with • = 0. 
Moreover, one readily verifies for a tridcndriform algebra (T, <, >, •) that (-Dt, >~), where -<,:= < + • 
and ;^:=>, is a dendriform algebra |Eb02j . 

Example 2. The summation operation ^ for 6* = 1 on a suitable algebra F of functions provides a natural 
example for such a tridcndriform algebra Tp = {F, <,>,•): 

{A < B){x) := A{x) ■ S{B){x), {A > B){x) := S{A){x) ■ B{x), {A • B){x) := A{x) ■ B{x). 
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One verifies without problems the tridcndriform relations. For instance, the first relation simply encodes: 

• S{B){x)^ ■ S{C){x) = A{x) ■ (5(5 • S{C)){x) + S{S{B) ■ C){x) + S{B ■ C)(a;)) 

identity ([8]). The associative product in this tridcndriform algebra then writes: 

{Ai.B){x) = A{x) ■ S{B){x) + S{A){x)-B{x)+A{x)-B{x). 

Within this example the reader may want to convince himself that Dtp = {F,^,,'^), where ^,:= < + •, 
i.e.: 

{A ^, B){x) = {A< B){x) + [A • B){x) 
and ;^:==>, is a dendriform algebra. 

One feels that dendriform algebras provide an elegant setting for a refined encoding of fundamental 
structures underlying integration and summation operations. In fact, further below we will show that this 
is true for a much larger class of operators, i.e. characterized by identity ([8|). 

Let A ^ A® k.\ be our dendriform algebra augmented by a unit 1: 

(22) a ^ 1 := a =: 1 )~ a l-(a:=0=:a^l, 

implying a^l ~ 1 * a ~ a. Note that 1*1 = 1, but that 1^1 and 1 >- 1 are not defined [RonOOj . [Cha02| . 
We recursively define the following set of elements of A[[A]] for a fixed a G A: 

w^^\a) :- 1 =: w'^^\a), w^J^\a) := a < {w^r^\a)), w^;^\a) := {w^y-'\a)) ^ a. 

Let us define the exponential and logarithm map in terms of the associative product (|16p . cxp*(a;) := 
Sn>o^*"/"''' 1°§*(-'- + ~ X]n>o(~-'^)"^*"/'^' respectively. In the following we first give a recursive 

expression for the logarithm of the solutions of the following two equations for a fixed a £ A: 

(23) X = 1 + Xa < X, Y y Xa. 

in A [[A]], in terms of the left pre-Lie product >. This will in particular encompass the Magnus expansion 
for the logarithm of a solution of a linear first-order homogeneous differential equation in a noncommutative 
algebra |Mag54| . Then we present the solutions of the two equations as an infinite product expansion of the 
exponential, which encompasses Fer's solution to a linear first-order homogeneous differential equation |F58j . 

Remark 1. For later use we may also write: 

(X - 1) = Aa + Aa ^ (X - 1) resp. {¥ ~ 1) = -Xa - (r - 1) ^ Xa. 

3. The pre-Lie Magnus expansion 
Formal solutions to ([23]) are given by: 

X = J2 A"wf7^(a) resp. Y ^ ^(-A)"w^"^(a). 

n>0 n>0 

Let us introduce the following operators in (A, -<,>-), where a is any element of A: 

L^[a]{b) := a<h [a] (6) := a 6 R^[a]{h) h ^ a i?^ [a] (6) := 6 a 

L<][a](6) := a<\h L^[a\{b) := a[> b R^[a]{h) := h <\ a Rt^[a\{b) b t> a 
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Theorem 2. Letfl' := Q'{Xa), a £ A, be the element of \A[[X]] such that X = exp*{Q') andY = exp*{-fl'), 
where X and Y are the solutions of the two equations \23\) . respectively. This element obeys the following 
recursive equation: 

(^"l - i-e,pff.|ni) '^°' = £'-'''"^"-'"'l''''^°'' 

or alternatively: 

where the Bi 's are the Bernoulli numbers. 

Proof. Let us notice that can be immediately derived from thanks to Lf^[b] = — i?<|[6] for any b G A. 
We prove p4|) . which can be rewritten as: 

(26) Xa= {n(Xa)). 

Given such fl' := il'{Xa) G AA[[A]] we must then prove that X := cxp*{n'{Xa)) is the solution of X = 
1 + Xa ~< X, where a is given by (|26p . Let us first remark that: 

(27) R^n'] ^ R^[n'] - Ly[n'], 

and that the two operators R^[n'] and i^[ri'] commute thanks to the dendriform axiom (fT4)l . We have 
then, using the three dendriform algebra axioms: 

= [ exp{sLy[n'])exp{-sR^[n']){n')ds 
Jo 

= [ exp*(srj') >- r2' ^ cxp* {-sil') ds. 
Jo 

Aa -< X = [ (exp*(sl7') y n' cxp*(-.sl7')) -< cxp*{n') ds 
Jo 

= [ cxp* {sil' ) y n' -<exp*{{l~ s)n')ds 
Jo 



So we get: 



(28) 



?i>0 p+q=n 

An iterated integration by parts shows that: 

p\q\ 



1 -\n\ 



{l-sYsP ds = 



(p + g+1)!' 



which yields: 



^ (n + 1)! ^ 

ri>0 p+(j=n 



On the other hand, we have 

(n + iy 



(29) X-l=exp*(»0-l = V , ,. »^*"+^ 

^ ^ ' ^ ^ ^ (n+ 1) 

n>0 

Equality (|26p follows then from the identity: 
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which is easily shown by induction on n. Analogously, one readily verifies that: 

y Aa = [ cxp*{-n') y {exp*{sn') y n' -<exp*{-sn'))ds 
Jo 

= / exp*((s - l)n') exp*{-sn') ds 

Jo 



= y y i-iYp+i^n'*'' yn' ^ / ds = y '^j^Krr y y^' 

^ ^ ' J p\q\ ^ (n+1)! ^ 



.... 

n>0p+q=n ""^ ^ n>0 ^ ' ' p+q=n 

□ 



Remark 3. It seems appropriate at this point to justify notation. We have chosen to write ft' to remind 
the reader of the fact that in the particular context of Example [T] (see subsection 15.31 below) one readily sees 
that n'{t) = n{t), see eq. 

Remark 4. In |EMP07bj we were able to show, using Hopf algebra and free Lie algebra techniques, i.e. the 
Dynkin idempotent map that: 

n'ia) = f (C{s) + Y^{-ir^ad^^){C{s))]ds, 
Jo V „>o / 

where C{t) = E„>o ^>H^"Ha)i"-'- 



4. The pre-Lie Fer expansion 

Let us come back to the dendriform power sums expansions (|23p: 

X = 1 + Xa-(X and Y = 1 - AY ^ a. 

We will mainly focus on the first one. Following Fer's original work |F58] we now make a simple Ansatz for 
its solution: 

X = exp*(Aa) *Vi 

and return this into the recursion. Recall the dendriform axiom (fT3|l . This then leads to the following: 

exp*(Aa)H<Vi = 1 + Aa ^ (exp*(Aa) * Vi) = 1 + A(a exp*(Aa)) ^ Vi. 

It is the goal to derive a recursion for Vi similar to the original one and then to iterate the process. 
Remember (1^^ and the dendriform axiom P^ : 

Vi = exp*(— Aa) + Aexp*(— Aa) * ^(a ^ exp*(Aa)) -< Vi^ 

= exp*(-Aa) + Aexp*(-Aa) -< (^[a -< exp*(Aa)) -< Vi^ + Aexp*(-Aa) >~ (^{a ^ exp*(Aa)) -< Vi^ 

= exp*(— Aa) + Aexp*(— Aa) -< ^(a exp*(Aa)) -< Vij + (exp*(— Aa) >- Aa ^ exp*(Aa)) -< Vi 

= 1 + (exp*(-Aa) - l) ^ (^1 + Aa -< (exp*(Aa) * Vi)'j + (exp*(-Aa) ^ Aa -< exp*(Aa)) ^ Vi 

= 1 + (^(exp*(-Aa) -1) < exp*(Aa)^ ^ Vi + (exp*(-Aa) >- Aa exp*(Aa)) -< Vi 

= 1 + (^(exp*(-Aa) - l) -< exp*(Aa) + exp*(-Aa) ^ Aa -< exp*(Aa)^ -< Vi. 

At this point we can repeat the above using the Ansatz: 

Vi :=exp*([/0*F2, 

where: 

(30) U[ := (exp*(-Aa) - l) exp*(Aa) + exp*(-Aa) >- Aa exp*(Aa). 
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In general we have: 



Vn :-cxp*(C/;j*K+i, 



with C/q := Aa and: 



(31) K (cxp*(-(7;_i) - 1) ^ cxp*([7;_i) + cxp*(-[/;_i) ^ (J^^i cxp*([/;_i). 
Such that we arrive at the following infinite product expansion for X: 

X = cxp*{Uo) * cxp*(C/{) * cxp*([/^) * ■ • ■ * exp*(t/;) * • • • . 
Analogously, using the Ansatz Y = Vi * exp*(— Aa) one shows that: 

r = • • • * cxp*(-C/;j * • • • * cxp*(-[/^) * cxp*(-C/() * cxp*(-[/^). 
Let us now examine a bit closer ()30p . By the foregoing calculation in Section [3] one readily verifies that: 

exp(i?<] [Aa])(Aa) = exp*(— Aa) :^ Aa ^ exp*(Aa). 
Lemma 5. Let (A, )^) be a dendriform algebra augmented by a unit 1 ^2^. Then, for a ^ A we have: 

(32) (exp*(— a) — l) -< exp*(a) = — / exp*(— sa) ;^ a -< exp*(sa) ds. 

Jo 

Proof. In the proof of Theorem [2] we have shown (see ([28]) . with —a replacing fl') that: 

exp*(— a) 1=1 cxp*(— sa) >- a ^ cxp*((s — l)a) ds. 
Jo 

This immediately yields: 

(exp*(— a) — l) ^ exp*(a) = — ^ / exp*(— sa) ^ a -< exp*((.s — l)a) ds^ -< exp*(a) 



exp*(— sa) ^ a -< exp*((s)a) ds 
by application of the first dendriform axiom. □ 
This lemma implies then for the general recursion ()3ip : 

(exp*(-;7;) - 1) ^ cxp*([/;) + cxp*(-c/;) ^ t/; ^ cxp*([/;) 
= cxp(i?^[L/;j)(c/,;)- / exp*(sc/;j ^ t/; -< cxp*(-s[/;)ds 



- ^^P{^<l[^n\){Un) n rrr, i (^/n) 

^ cxp(-L,ra)(t/;) + -P^-;->[^^])-\ ^;) 

(L^[£/;] + l)exp(-L^ [£/;])-! 

and one shows that this then gives the nice identity essentially encoding Fer's classical expansion in terms 
of a prc-Lic product: 

(exp*(-f/;) - 1) ^ exp*([/;) + exp*(-C/;) exp*([/;) = ^ 

i>o ^ 

In the following theorem we summarize the above by formulating Fer's expansion ([5|). 
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Theorem 6. Let {A, )^) be a dendriform algebra augmented by a unit 1 i22\) . Let Uq :— Xa, :~ U'^{a), 
n G N, a G A, be elements in \A[[X\\, such that 

X^\{* exp*(f/;) Y = n* cxp* (-[/;) 

Ti>0 n>0 

where X and Y are the solutions of the equations \2S\) . Then these elements obey the following recursive 
equation: 

(33) t^;+i:=E7^TlT!^'>[^"]^'^^^"^' 

The presentation given here in the context of dendriform algebras reduces to Fer's classical expansion 
when working in the dendriform algebra of Example [T] Let us recall that Iserles |Ise84| rediscovered Fer's 
result calling it the method of iterated commutators. Munthe-Kaas and Zanna |MZ97| further developed 
Iserles' work in the context of Lie group integrators. 

5. Applications 

In this section we collect some applications where the foregoing straightforwardly implies known as well 
as new results. 

5.1. Associative algebras. Any associative algebra {A, *) can be seen as a dendriform algebra with -<= * 
and ^=0 (or alternatively -<= and y= *). In this case the pre-Lie operation < reduces to the associative 
product, and equation reduces to: 

(34) \a = l-cxp*{-n'), 

hence il' = — log*(l — Xa). Its exponential X = {1 — Xa)*^^ = 1 + Xa + X^a * a + ■ ■ ■ indeed verifies 
A = l + Aa*A. 

5.2. Rota— Baxter algebras. Recall jBax60| IEb02| IRot69| that an associative Rota-Baxter algebra (over 
a field k) is an associative fc-algebra A endowed with a /c-linear map R : A ^ A subject to the following 
relation: 

(35) R{a)R{b) = R{R{a)b + aR{b) + Bab) . 

where 9 ^ k. The map R is called a Rota-Baxter operator of weight 9. The map R := —Oid — R also is a 
weight Rota-Baxter map. Both the image of R and R form subalgebras in A. Associative Rota-Baxter 
algebras arise in many mathematical contexts, e.g. in integral and finite differences calculus, but also in 
perturbative renormalization in quantum field theory jEGM06] . 

A few examples are in order. On a suitable class of functions, we define the following Riemann summation 
operators 

[x/0\ [x/e]-i 

(36) Re{f){x) := Bf{n9) and R'g{f){x) := ^/("^)- 

n—1 n—1 

[x/e\ [x/e\ [x/e] [x/e] [x/e\ 

Y9f{nB)\r£eg{ni9)\^i E + E + E )^'/(n%(m0) 

n—1 ■'^ ^ ?n — 1 ^ 71 > ni — 1 ni > 7i — 1 ni — n — 1 ■^^ 

[x/O] , m s [x/e] n s [x/e] 

= E eHYf{k9)\g{m9) + E eMYg{k9)\f{n9) - E e^f{ne)g{n9) 

m=l ^fc=l ' n=l ^fe=l ^ n=l 

(37) = Re{Reif)g)ix) + Reif Re{g)){x) + 9Re{fg){x). 

Similarly for the map R'q. Hence, the Riemann summation maps Rg and R'g satisfy the weight —0 and the 
weight 9 Rota-Baxter relation, respectively. 



Observe that 
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Let US give another example, very different from summation and integration maps. Let A be a K- 
algebra which decomposes directly into subalgcbras Ai and A2 , A ~ Ai (B A2, then the projection to Ai, 
R: A A, R{ai, 02) ~ ai, is an idcmpotcnt Rota-Baxter operator, i.e. of weight 6 = —1. Let us verify this 
for a, 6 e A = Ai ® A2 

R{a)b + aR{b) - ab = R{a){R{b) + {id - R){b)) - {R{a) + {id - R){a)) {id - R){b) 

= R{a)R{b) - {id- R){a){id- R){b) 

such that applying R on both sides kills the term {id — R){a){id — R){b) without changing the term R{a)R{b), 
as R{id — R){a) = since Ai, A2 are subalgebras. 

Proposition 7. |Eb02| Any associative Rota-Baxter algebra gives a tridendriform algebra, (Tr, <, >, •g), 
in the sense that the Rota-Baxter structure yields three binary operations: 

a < b := ai?(&), a > b := R{a)b, , a •gb := 6ab, 

satisfying the tridendriform algebra axioms i2S^) . 

The associated associative product *g is given by 

a *e 6 aR{b) + R{a)b + 0ab 

It is sometimes called the "double Rota-Baxter product" , and verifies: 

(38) R{a*gb) ^ R{a)R{b), R{a *e b) ^ -R{a)R{b) 

which is just a reformulation of the Rota-Baxter relation ([35]) . It follows from the general link between 
dendriform and tridendriform algebras that any Rota-Baxter algebra gives rise to a dendriform algebra 

structure, {Dfj,^,"^), given by: 

(39) a-<b := aR{b) + Oab = -aR{b), a >- b := R{a)b. 
For completeness we mention the following. 

Remark 8. Rota-Baxter Dendriform algebra It is easy to verify that {Dp, y-) defines a Rota-Baxter 
dendriform algebra with weight 6 Rota-Baxter map R : Dr Dr, that is: 

R{a)-<R{b) = R{R{a) b + a R{b) + Oa b) 
R{a)>-R{b) = R{R{a) ^ b + a ^ R{b) + Oa >- b) 
This provides an example for a non- associative Rota-Baxter algebra. 
The (weight 9) Rota-Baxter (left) pre-Lie operation corresponding to (|18p is given by: 

(40) a>b = R{a)b~bR{a)-eba^[R{a),b]-eba. 

Let us remark that the underlying vector space A equipped with associative product *g is again a Rota- 
Baxter algebra with weight 9 Rota-Baxter map R. Whereas, the pro-Lie algebra {A, [>) is a Rota-Baxter 
pre-Lie algebra with weight 6 Rota-Baxter map R, i.e. another example for a non-associative Rota-Baxter 
algebra. For completeness, recall the known fact that, analogously to ([55]) . we can define a new pre-Lie 
product [>£/ on {A, o), and {A, \>g) is again a Rota-Baxter pre-Lie algebra with weight 9 Rota-Baxter map 
R. 

Notice that if we suppose the algebra A to be unital, the unit (which we denote by 1) has nothing to do 
with the artificially added unit 1 of the underlying dendriform algebra. We extend the Rota-Baxter algebra 
structure to A by setting: 

(41) R{1) := 1, ^(1) := -1 and 1.x = x.l = for any x £ A. 

This is consistent with the axioms ()22p which in particular yield 1 >- x = R{l)x and x ~< 1 ~ —xR{l), in 
coherence with ([39]) . 
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Now let US suppose that the Rota-Baxter algebra A is unital, and introduce the weight 9 G k pre-Lie 
Magnus type recursion, fig Q'g{Xa) G AA[[A]], where the Rota-Baxter operator R is naturally extended 
to A[[X]] by fc[[A]]-hnearity: 

(42) %{\a) = J2^LM]^"^H^a), 

m>0 

Theorem [5] together with relation ([55]) implies for a fixed a G A the following corollary. 

Corollary 9. The elements X := -R{X) = exp{-R{n'g{Xa))) and Y := R{Y) = cxp{-R{n'g{Xa))) in 
A[[X]] solve the equations: 

(43) X ^1- XR{aX) resp. Y = 1- XR{Ya). 

Proof. Recall the link between dendriform and tridcndriform algebras, and the Rota-Baxter relation, i.e. 
relations ([7]). Using p8|) . we have: 

X = -R{X) 

= -R{1 + Xa-<X)^ -R{1 - XaR{X)) = -R{l + XaX) 
= l-XR{aX). 

and similarly: 

Y = R{Y) 

= R{1 -XY )-a) = R{1 - XR{Y)a) = R{1 - XYa) 
= 1 - XR{Ya). 

□ 

We may summarize the picture we developed so far in the following diagram relating the pre-Lie Magnus 
expansion to the original Magnus expansion as well as to the Spitzer's identity: 



\ ^m>0 'J 




exp(7?(a)) 



^—0, com. 

On the left wing we regard the limit d zero leading to Magnus' expansion, whereas on the left wing 
we observe the reduction of the pre-Lie Magnus expansion to the logarithm in Spitzer's identity when the 
underlying algebra is commutative. Finally both formulas reduce to the classical exponential solution in a 
commutative weight zero Rota-Baxter algebra. We refer the reader to [EGM06j for more details. 

We define here the weight 6 G k pre-Lie Fer type recursion, 1/,^ g := 1/,^ g{a) € A^[[A]]: 
(44) K+iA^a) - J2t^LM,of \Ke)^ " > «■ 

By the same line of arguments Theorem [6] implies the next corollary. 
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Corollary 10. The element 

X := n cxp(-i?(C/;,(a))) Y := [f exp(-i?(C/;,(a)) 



ri>0 n>0 

in ^[[A]] solves the above recursions. 

Associative Rota-Baxter algebras are essentially characterized by a natural factorization theorem related 
to the equations (|43)) . see Atkinson [At63j for more details. 

Recall |EGM06j that there is a unique (usually non-linear) biiection xe '■ AA[[A]] AA[[A]], which is a 
deformation of the identity, such that for any a G AA[[A]] we have 3: 

exp(— 6'a) = exp(i? o xe{ct)) exp(i? o xe(«)) • 

The so-called BCH-recursion map xo is recursively given by: 

Xe{a)^a + ^ BCH(i? oxe{a),Ro xe{a)) , 

or alternatively: 

= a + i BCH(^?a, R o xeia)) , 

where BCH(a;, y) is the Baker-Campbell-Hausdorff series defined by exp(a;) cxp{y) — exp(x+y + BCH(x, y)) . 
Recall Atkinson's theorem jAt63| : 

Y{1 ~e\a)X = 1, 

hence Y^^X^^ = 1 — dXa. We deduce immediately from the very definitions of X and Y that we have: 

(45) 1 - 6lAa = exp(-6lae) = ciip{R{n'g)) exp(^(f7^)) , 
with ag := ag{Xa) :~ — ^ log(l — 6Xa). We then infer from ((45|) the equality: 

(46) % . %iXa) . xeiae) = Xe {-'^^^^^j^) , 

which was conjectured in |CEFG06] in a similar context. From we get for any a E XA[[X]]: 

(47) Xe{a)^n',{l^^^^- 

5.3. The weight zero case and the classical Magnus expansion. For a weight 9 — 0, Rota-Baxter 
algebra the pre-Lie product pO|) reduces to: 



(48) a[>b = [R{a),b]. 

This simplifies the weighted pre-Lie Magnus type recursion (|42| to the classical Magnus recursion: 



n'„{Xa) ^ XoiXa) = ^ :^L^[;^o](™)(Ac 



m! 

m>Q 



(49) - E^-d^nLi^-y 

m>0 



An important example for a weight zero Rota-Baxter algebra is any algebra F of operator- valued functions 
on the real line, closed under integrals , say smooth n x n matrix- valued functions. Recall example ([T]), 
showing that is a dendriform algebra. The associative product then writes: 

X X 

A * B{x) := A{x) ■ [ B{y) dy + [ A(y) dy ■ B{x) 



^Notice the sign change compared to |EGM06] . This is due to a conventional sign change in the definition of the weight: 
an idempotent Rota-Baxter operator is now of weight —1 with the new convention. The two recursive formulae for Xfl which 
follow thus differ a little bit from those given in the aforementioned reference. 
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such that: ^ ^ 

exp(^ A{y)dv)^ (cxp*(A(y))) dy. 

Recall our convention (|¥T|) for the dendriform unit. Then the classical Magnus recursion ([2]), ^{t) = Vl{A){t), 
ioi R:— reads: 

m>0 m>0 

Similarly, in the case of the weighted pre-Lie Fer expansion as expected we recover the original formula. 

We showed in reference jEGM06j that the recursion xe reduces to the classical Magnus recursion p9)) in 
the hmit 6* ^ 0: 

Xe(Aa) Xo(Aa). 

We recover this result from (j47|) . ((42|) . and from the fact that ag{\a) Aa when 6 ^ Q. Hence, we have 
proven that the weight 9 BCH-recursion xe derives naturally from the classical Magnus recursion via the 
non-linear change of variable ag. 

Let us finish this brief note with an interesting observation to be further explored in a future work. In 
the context of Example [1] and the initial value problem ^ we have seen that the classical Magnus and Fer 
recursion can be rewritten using the (weight zero) pre-Lie product introduced by the integral operator j^: 

= E ^L^mA){s)r^\A{s)), 

m>0 

and 

Un+i.o{A){s) =Y,j-^^L^[U,,,,{A){s)]^'\UnAA{s))). 

l>0 ^ ^' 

It seems to be natural to ask whether the extra pre-Lie structure in the case of the Magnus as well as 
the Fer expansion implies a possible reduction in the number of terms. Here we would like to indicate that 
this is the case, by explicit verification up to fifth order in the case of the Magnus expansion, using the 
dendriform algebra respectively the induced pre-Lie product presented in Example [TJ We should emphasize 
the fact, that the following applies, of course, to the more general weight 6 pre-Lie Magnus and Fer type 
recursions l^^ . ([33]), respectively. 

We introduce a dummy parameter A for convenience. Obviously, as the (left) pre-Lie relation p8|) is a 
ternary one we expect it to be available only from third order upwards in this parameter. 

(50) n{t)= f nQ{AX){s) ds = \ f A{s)ds-\^l f A[>A{s)ds 
Jo Jo 2 Jq 

+ A' (A o (A o A)) (s) + i ((A > A) > A) (s)) ds 

+ (^-l{iA:>A):>A)> A{s) - 1 (A > (A > A)) > A(s) 

- ^(^Ai> {{At> A)t> A){s) + iA\> A)\> {At> A){s)^^ ds + 0(5) 

Recall that A o A{s) = [J^ A{u)du, A{s)]. We see that at third order no further reduction of terms is 
possible. At fourth order we find a reduction to two terms using the pre-Lie relation (jlSp . Indeed, one 
verifies that: 

^{{Ao A) i> A) t> A+ ^(^{Ai> {Ai> A)) > A + A> {(A^ A) c> A) + {A\> A) \> (A o A)) 
using that, thanks to the pre-Lie relation: 

{Ac> A) ^ {A[> A) ^ {{A> A) > A) t> A - {A> {A[> A)) > A + A\> {{A[> A) [> A) 
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equals: 

1({A[> A)[> A) \> A + ^A\> ({A \>A)\>A). 
6 12 

At fifth order we observe a reduction in the number of terms from ten to seven. More details and a complete 
analysis of this apparently new structure in the Magnus (and Fer) expansion will be provided in a forthcoming 
work. 

Following the seminal work of Iserles and N0rsett [IN99| , using planar rooted binary trees to encode the 
combinatorial setting in Magnus expansion, we may present the foregoing calculation more transparently. 
The binary tree: 

- A[>A. 

At fourth order we have: 




which reduces to: 




thanks to the left prc-Lic relation: 
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